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Page  iv 

paragraph  3,  second  line:- 
replace  "minimum"  by  "maximum" 

Page  6 

equation  (8)  shall  be  written  as: 

%i-f§ 

Figure  9 

graph  No.  41  (  >10*1.63)  shall  be  replaced 
by  graph  No.  41A  (March  5,  1963) 

Page  35 

paragraph  2,  second  line:- 
Figure  12  shall  read  Figure  13 

Two  seta  of  equations  in  Part  II  have  been  numbered  15  to  22, 
one  set  on  pages  21  to  23  and  the  other  set  on  pages  24  to  27. 


March  6,  1963 
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SUMMARY 


The  temperature  difference  between  the  bright  and  dark  aide  of  the 
apecified  800  -  ft. ,  half  inch  diameter,  Beryllium  Copper  element, 
when  exposed  to  sun  radiation,  will  reach  12°F.  This  difference  will 
cause  the  element  to  bend  away  from  the  radiation  source  in  an  asymptotic 
form.  The  tip  deflection  will  be  about  550  ft.  ,  with  a  minimum  radius  of 
curvature  of  333  ft. 

The  location  of  the  point  of  minimum  radius  of  curvature  can  be 
anywhere  along  the  axis  of  the  element  from  the  root  to  the  tip,  depending 
on  the  element's  initial  direction  with  respect  to  the  radiation  source. 

The  point  of  minimum  radius  of  curvature  coincides  with  the  point 
of  minimum  element  temperature ,  which  will  be  about  850°R  with  sun  and 
earth  radiation,  and  785°  R  in  outer  space. 

The  temperature  difference  along  the  cross  section,  and  the  re¬ 
sulting  profile  change,  are  determined  by  the  material  and  its  surface 
absorptivity  only  and  are  independent  of  the  emissivity;  while  the  maximum 
temperature  depends  on  the  ratio  of  absorptivity  to  emissivity  only. 

The  magnitude  of  the  time  constant  of  the  specified  element  with 
respect  to  step  changes  in  radiation  is  about  3  minutes.  The  final  profile 
(96%)  will  be  established  in  about  7  minutes  after  a  sudden  extension  of 
the  element  into  the  radiation  field. 
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EVALUATION  AND  DISCUSSION 
OF  ANALYTICAL  RESULTS 

1.0  INTRODUCTION 

! 

This  study  of  the  effect  of  rediant  energy  impinging  on  long  tubular 
elements  in  space  investigates  the  temperature  distribution  along  the 
perimeter  and  the  axis  of  the  element  and  the  resulting  profile  changes. 

The  results  are  condensed  into  equations  and  into  a  Profile  Chart,  allowing 
examination  of  the  particular  profile  of  an  element  in  outer  space ,  or  in 
an  orbit  around  the  earth.  The  time  constant  of  the  element  is  also  in¬ 
vestigated  analytically. 

1 . 1  Specification  of  the  Element 


The  specified  element  is  considered  as  a  closed  tube  for  most  of  the 
investigations.  A  correction  factor  is  discussed  and  proposed  for  an  open 
tube  with  180°  overlap.  Numerical  values  of  temperatures  along  the  per¬ 
imeter  and  the  axis ,  and  values  for  the  minimum  radii  of  curvature  and 
tip  deflections  are  calculated  for  an  element  with  the  following  specifications: 


Element  Material 

Element  Diameter 

Element  Length 

Wall  Thickness 

Material  Density 

Thermal  Conductivity 

Thermal  Expansion  Coefficient 

Absorptivity  Ratio 

Emissivity  Ratio 

Specific  Heat 

Overlap  Factor 


Beryllium  copper 
d  =  0.  04167  feet  (0.  50  inches) 

1  =  800  feet 

t  =  0.000167  feet  (0.002  inches) 
/>  =  500  lb /ft 3 
k  =  44  BTU/ft  -  °F  -  hr 
e  =  1 . 04  x  10*3 
<X  =  .45 
€  =  .  10 

C  =  .10  BTU/lb  -  °F 
?  =  1.5 


1  •  2  Specification  of  the  Radiant  Energy 


Sun  and  earth  radiation,  and  reflection  of  sun  radiation  by  the  earth 
and  the  atmosphere,  are  taken  into  account.  These  are  assumed  to  have 
the  following  intensities  :  - 
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Sun  Radiation 
Earth  Radiation 
Earth  Reflection 


460  I1T(T/HR 


FT 


2 


70  DTI! /HR  FT2 
160  BTU/HR  FT2 


Sun  radiation  ie  coneidered  to  be  constant  in  the  scope  of  this  analy¬ 
sis.  Earth  radiation  and  reflection  are  subject  to  correction  for  altitude, 
the  values  applied ,  referring  to  the  upper  layers  of  the  mesosphere . 
which  ie  assumed  to  be  500  miles  above  the  surface  of  the  earth  . 


The  amount  of  heat  entering  the  element  is  also  a  function  of  the 
incident  angle.  This  angle  is  defined  in  the  analysis  as  the  angle  between 
the  radiation  direction  and  the  axis  of  the  element  as  seen  from  the  root 
of  the  element. (See  diagram  below.) 


2.0  THE  ELEMENT  IN  OUTER  SPACE 

2. 1  The  Maximum  Temperature  Difference  Along  the  Element  Penmoter^ 

The  maximum  temperature  difference  between  the  bright  and  dark 
side  of  the  element  occurs  when  the  incident  angle  90  .  It  is  obtained 

by  the  equation: 

AT»*i?irCXs',sSmes  <l> 

It  is  independent  of  the  emmissivity  ratio  £.TS  (max>  »8  lor  the 

specified  element. 
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2-2  The  Temperature  Profile  Along  the  Element  Perimeter 

The  Temperature  Profile  is  shown  in  Figure  2.  It  is  symmetrical  on 
both  sides  of  the  diameter  from  the  bright  side  to  the  dark.  It  is  obtained 
from  the  equations: 

A  Ts  M  -  ATs  m,x  [  f+  «#  JL  -  (2, 


a  (1;  *  ats  ^a„  -f-)1-  -ssj 


Equation  2  gives  the  profile  in  the  first  and  fourth  quadrants,  and  equation 
3  the  profile  in  the  second  and  third  quadrants.  In  both  cases,  X  is  a  point 
on  the  perimeter  expressed  in  terms  of  IT  X  IT 

2.3  Element  Profile  Change 

The  temperature  difference  between  the  bright  and  dark  sides  of  the 
element  will  cause  the  element  to  bend  away  from  the  radiation  source. 

The  profile  of  the  bent  element,  the  tip  deflection,  and  the  radii  of  cur¬ 
vature  along  the  element  axis,  are  determined  by  a  "Profile  Number" 
which  is  defined  as 


Profile  Number 
It  has  a  value  of  3  x  10'^  for 


the  specified  element. 


Using  this  profile  number,  the  element  shape  can  be  derived  by  the 
equation: 


tr>  cos  (Y*N») 


This  profile  has  been  plotted  in  Figure  2  for  the  specified  element 
and  a  radiation  incident  angle  ©s  =  90°.  The  X  and  Y  axes  are  in  feet; 
and  the  length  of  the  element  has  been  traced  on  the  plotted  profile. 
Figure  3  shows  the  same  element  for  a  radiation  incident  angle  ©5  =  135° 
Figure  4  shows  the  profiles  of  a  dipole  element  with  incident  angles  from 
0  to  1 80° . 


Explanation  of  the  Profile  Chart 


The  profile  curves  of  Figures  2,  3  and  4  are  extracted  from  the 
ProfileChart  shown  in  Figure  5.  Any  800  foot  length  traced  on  the  curve 
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with  the  Profile  Number  3  x  10*3,  can  be  a  profile  of  the  specified  element. 
The  root  of  the  element  is  located  at  the  intersection  point  of  the  incident 
angle  line  with  the  profile  curve,  and  the  tip  of  the  element  is  found  by 
tracing  the  length  of  the  element  (800  feet)  clockwise  on  the  profile  curve. 

The  Profile'Chart  and  equation  4  for  the  Profile  Number,  imply  that 
elements  with  different  design  parameters  can  have  the  same  profile  change 
as  long  as  the  design  parameters  in  equation  4  yield  the  same  Profile 
Number.  Hence  one  can  examine  the  change  of  profile  for  a  change  of  a 
particular  design  parameter,  by  comparing  the  profile  curves  with  the 
appropriate  number s. 

2 . 5  Several  Elements  on  One  Satellite 

Ao  long  as  the  design  ensures  that  the  incident  radiation  angle  is  the 
sam-,  for  all  the  elements,  several  elements  on  one  satellite  will  have  the 
same  profile.  A  crossed  dipole  of  elements  as  specified,  will  have  a  con¬ 
figuration  as  shown  in  Figure  6.  The  configuration  will  be  maintained  in 
the  orbit  around  the  sun  as  shown  in  Figure  7.  The  elements  have  the 
tendency  to  align  themselves  and  the  satellite  radially  with  respect  to  the 
orbit,  so  that  the  same  side  of  the  satellite  surface  always  faces  the  Sun, 

The  backward-bent  elements  which  are  exposed  to  the  sun's  pressure,  will 
create  a  correcting  torque  as  soon  as  the  satellite  deviates  from  the  equilib¬ 
rium  either  in  the  pitch  or  yaw  axis.  The  dynamic  properties  of  the  elements; 
their  damping  effects,  and  time  lags,  have  to  be  included  in  the  orientation 
and  stabilization  system  of  the  satellite. 

2 . 6  Tip  Deflection 

The  value  X  in  Equation  5  represents  the  tip  deflection  from  the 
initial  straight  profile.  It  can  be  obtained  only  by  tracing  the  profile. 

However,  for  very  small  Profile  Numbers,  or  for  short  elements,  Y  can 
be  replaced  by  the  element  length;  £  and  Equation  5  can  be  simplified  to: 

(6) 

Thus  the  deflection  for  a  100  foot  length  would  be  15  ft.  ,  when  ex¬ 
posed  to  s»m radiation  only. 

2. 7  Minimum  Radius  of  Curvature 

The  greater  the  temperature  difference  on  the  cross  section,  the 
more  the  element  will  bend  at  that  section,  and  the  smaller  will  be  the 
radius  of  the  curvature.  The  location  of  the  minimum  radius  of  curvature 
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will  be  at  the  point  where  the  incident  angle  is  90°.  In  Figure  2  the 
location  is  at  the  root  of  the  element  and  in  Figure  3  wheL  the  tangent  on 

U  “  ,h'  Th.  ,«ilu.  ..  any* point 

Rf*> =  TT  sin  (?) 

where  0ft) is  the  angle  between  the  tangent  and  the  initial  direction.  The 
minimum  radius  of  curvature  for  the  specified  element  is  333  feet,  .when 

comnut  f, ,Un  “dlation-  An«le  ®^)can  be  obtained  from  the  graph  or 
computed  from  the  equation:  *  ** 


2  •  8  Maximum  Element  Temperature 

The  element  temperature  depends  on  the  absorptivity  arid  emissivitv 

,„cW.„  „dl.,lon 

:,::irz::;rrp/uTr:;yA'  *w  p“tat  x  *'“«  -•  *-•  -•  — ... 


T(») 5  (l&W  w3*sin  e<*>)* 


where  6u)is  again  the  angle  between  the  tangent  on  the  element  and  the 
radiation  direction.  The  location  of  the  point  of  maximum  temperature 

Figure  Y  h  l04Cati0n  °f  the  P°int  of  minimum  radius  of  curvature 

and  Fia  8lTY  ma5Un?Um  element  temperature  versus  the  ratio  <*-/e 
and  Figure  9  depicts  the  temperature  along  the  axis  of  one  element  fn  two 
different  positions  with  respect  to  the  sun.  The  maximum  temperature  of 
the  specified  element  is  785°R  when  exposed  to  sun  radiation  only 

2 . 9  Element  Response  Time 

It  has  been  assumed  that  the  element  will  have  found  structural  and 

,he™*‘  wh'»  *»•  >•«*»  of  hi,,,..,  t.mp„r,:rr.» 

the  tim.'i  i  ‘  condition-  the  ‘ime  constant  Ct  .  which  is  defined  a. 

after  t0  reach  63-2%  of  maximum  temperature 

rise  after  a  small  step  increase  in  radiant  energy,  is:  - 
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7  c  'Smq* 

"C  "  for  the  specified  element  would  be  3  minutes. 

However,  the  actual  temperature  response  curve  for  a  large  step 
increase  involving  sudden  extension  of  the  element  into  the  sunlight,  must 
be  derived  from  the  equation: 


^TfstepJ  *  -^§p(Ts\„a*“T^)  ATT  (step) 


The  transient  temperature  profile  is  shown  for  the  specified  element 
in  Figure  10.  The  element  needs  7  minutes  to  reach  96%  of  the  final  profile. 


3.0  THE  ELEMENT  EXPOSED  TO  RADIANT  ENERGY  FROM  SUN  AND 

EARTH 

3-1  Analysis  and  Definition  of  Earth  Radiant  Energy 

A  satellite  in  orbit  around  the  Earth  will  be  equally  exposed  to  earth 
radiation  due  to  the  earth  temperature  at  all  positions  of  its  orbit.  The 
earth  radiation  constant  has  to  be  corrected  for  the  distance  of  the  satellite 
from  the  earth. 

Reflected  radiant  energy  can  reach  the  satellite  only  as  long  as  the 
sun  and  earth  are  in  opposition  with  respect  to  the  satellite.  However, 
since  the  reflected  sunlight  is  diffused  in  the  atmosphere,  it  has  been 
assumed  that  the  reflection  approaches  zero,  when  the  angle 7£  between 
sun,  Earth  and  Satellite  reaches  120°.  (See  diagram  below.)  The  variation 
of  reflected  radiation  with  7}  ,  between  0  and  120  degrees,  can  be  repre¬ 

sented  by  a  complicated  expression  derived  in  Reference  1.  However,  it 
has  been  found  that  the  variation  can  be  reasonably  approximated  by  the 
function  cos  (.75  7^).  a. 

* — 


i 


Both  radiant  energy  sources  -  earth  radiation  and  earth  reflection 
are  subject  to  correction  for  altitude  with  the  factor 


since  radiation  intensity  varies  inversely  with  the  square  of  the  distance. 
R  is  the  radius  of  the  earth,  and  A  is  the  altitude  of  the  satellite  above 
the  earth  surface. 

The  effect  of  radiant  energy  from  the  earth,  also  depends  on  the 
magnitude  of  the  angle  between  the  sun  and  earth  radiation  with  respect 
to  the  element  cross  section.  This  angle  is  illustrated  in  the  diagram 
below  and  should  not  be  confused  with  the  angle  ^  .although  both  may 
have  the  same  value  in  particular  positions. 


3.2  The  Maximum  Temperature  Difference  Along  the  Perimeter 

The  maximum  temperature  difference  and  the  profile  along  the 
perimeter  is  obtained  by  multiplying  the  appropriate  equations  1  to  3  by 
the  factor  B  to  get: 


'SRE  4tk 


as3ss'mes*B 


wherein  B  represents  the  effect  of  earth  reflection  of  sunlight  and  earth 
radiation.  B  is  developed  as  follows:  the  total  amount  of  impringing  radia¬ 
tion  which  will  cause  the  temperature  rise  consists  of: 

Sun  radiation:  a3sSiri0s  lc 

Earth  reflection:  cOrCOS  .IS'H*  Sin0t*  ra  COS  © 

Earth  radiation:  £  Jr  Sin  ©,  x  Fa  COS*  S 

The  angles  9S  and  0fcwere  defined  in  1.2  and  the  angle  .7511  as  well  as 

the  factor  F  in  3 . 1 .  Thus  the  total  amount  of  impringing  radiation  is:  - 
A 


J 
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I3sa3ssines  ±|cORcos  .157)+  e3e)  sin9E  xFa  cos lS  <14> 
I 3  =oQs  Sine^it (-2&-  cos  .157}+  -l£—  FAcos*Sj  ( 1 5) 

where  the  sign  +  applies  for  6<  ^  and  the  sign  —  for 
The  term  in  brackets  is  called  B  thus: 

B,  =  I±(.3S5C0S  .lS7}  +  .l55-^)fft§*FACOSl6  (16) 

and  when  reflection  approaches  aero  with  O'}  >  120° 

Faco5*6  (17) 


Inspection  of  the  terms  for  B  reveals  that  the  magnitude  of  B  depends 
mainly  on  the  magnitude  of  the  four  angles  which  define  the  position  of  the 
element  with  respect  to  the  sun  and  earth.  This  is  because  tho  ratio 
is  very  small,  and  is  constant  for  the  specified  element. 


The  extreme  values  for  B  occur  in  the  positions  in  the  diagram 
below.  The  orbit  being  assumed  to  be  500  miles  above  the  earth’s  surface. 


es=90r 
©6»  u>5° 
7}  =  105° 

S~  0" 


es=ee-9o# 
71 « o° 
h  =  150* 
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3-3  The  Temperature  Profile  Along  the  Perimeter 

The  temperature  profiles  for  the  above  defined  positions  are  shown 
in  Figure  11.  They  represent  the  borderlines  for  all  possible  conditions. 
They  are  computed  from  the  equation: 

ATser=ATsx&  (js) 


Since  Bmjn  will  increase,  and  Bmax  will  decrease  with  orbits  at 
higher  altitudes,  both  will  approach  Unit  1  at  an  altitude  of  about  20,000 
miles.  This  means  that  the  influence  of  earth  radiation  and  reflection  will 
become  virtually  negligible  at  this  altitude  and  above.  The  temperature 
profile  then  will  be  that  indicated  by  the  dotted  line.  This  is  the  same  as 
shown  in  Fig.  1. 

3.4  Effect  on  Earth  Radiant  Energy  on  the  Element  Profile 

The  Profile  Number  of  the  element  will  be  changed  in  the  same  ratio 
as  the  temperature  gradient  along  the  perimeter  is  changed.  Thus  equation 
4  becomes: 


Nser=Ns*B  (19) 

Because  Bmilj  was  found  to  be  .  62  for  the  specified  element  with 
=  .22  the  minimum  Profile  Number  for  the  specified  element  when 
exposed  to  sun  and  earth  radiant  energy  will  be: 


Nser  (min) B  l-SfcxlO 


The  appropriate  profile  for  this  Profile  Number  is  marked  in 
Figure  12  on  the  Profile  Chart,  for  comparison  with  the  profile  for  the 
number  1  x  10"3.  it  is  the  lowest  Profile  Number  which  can  occur  in  the 
most  favourable  position  with  sun  and  earth  in  opposition,  with  respect  to 
the  satellite.  For  all  other  positions  particularly  for  all  higher  altitudes 
the  Profile  Number  will  increase  until  it  approaches  the  value  3  x  10'3  at 
an  altitude  of  about  20,000  ft. 

3 . 5  The  Elements  in  Orbits  Around  the  Earth 

Unless  the  satellite  travels  across  the  shadow  of  the  earth,  where 
the  elements  will  have  almost  a  straight  profile  -  the  elements  have  the 
tendency  to  align  themselves  to  the  direction  of  the  sun  radiation.  This 
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applies  for  all  orbits,  polar  orbits  as  well  aB  for  equatorial  or  semi- 
equatorial  orbits. 

In  a  polar  orbit  (Figure  13) ,  with  the  satellite  orientated  to  the 
earth  so  that  always  the  same  surface  will  face  the  earth,  the  elements 
will  keep  their  configuration  with  respect  to  the  earth  and  to  the  sun. 

They  will  not  influence  the  orientation  and  stabilizing  system  of  the 
satellite. 

In  a  semi -equatorial  orbit  however  as  shown  in  Figure  14  the  elements 
will  try  to  align  themselves  and  the  satellite  to  the  sun  radiation,  so  the 
satellite  will  perform  one  revolution  around  its  axis  while  performing 
one  orbit  cycle.  Thus  when  the  satellite  would  have  an  orientation  control 
to  face  the  earth  always  with  the  same  side  of  the  surface  the  orientation 
control  has  to  cater  for  the  disturbance  caused  by  the  elements. 

In  such  cases  the  elements  would  take  profiles  as  shown  in  Figure 
15  and  16.  In  Figure  15  only  one  element  is  assumed,  pointing  to  the 
earth  in  its  initial  position.  For  Figure  16  a  dipole  of  two  elements  is 
assumed  in  the  plane  of  the  orbit,  as  tangents  on  the  orbit  path  in  their 
initial  position.  This  configuration  assumes  that  the  orientation  control 
keeps  the  satellite  always  in  the  position  facing  the  earth  with  the  same 
side  of  the  surface. 

3. 6  Maximum  Element  Temperature  Along  the  Element  Axis 

The  temperature  of  the  element  on  the  point  of  Minimum  Radius  of 
Curvature  is  a  function  of  the  total  radiant  energy  impinging  on  the 
element.  Thus 


0so(  sin  9S 


in  equation  9  is  replaced  by; 


sin  005  «•"  ©s  c0& 


which  yields 

Tser  =TS  *  B* 


(20) 


but  with  always  the  sigh  i  after  1.  (Equation  16  in  3.2.) 
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Thus  the  maximum  temperature  occurs  for 
Fa  =  I  and  es«eE»90°  and  11-0*  and  S=lS0°, 

when  sun  and  earth  are  in  opposition  with  respect  to  the  satellite. 

In  this  case  the  maximum  temperature  for  the  specified  element  is: 

Tse#-(max.)=aS°*R 

3. 7  Piement  Response  Time 

Replacing  Tg  max  in  equation  10  by  Tser  =  Tfi  max  B’3/4  yields 
the  response  time  of  the  element  when  exposed  to  sun  and  earth  radiant 
energy  simultaneously.  It  would  be  about  2.35  minutes  instead  of  about 
3  minutes  so  that  the  overall  pattern  of  the  transient  curve  in  this  case 
will  be  similar  to  the  curve  as  shown  in  Figure  10. 


4.0  DISCUSSIONS  OF  BASIC  ASSUMPTIONS 
4-  1  Average  Temperature  and  Internal  Radiation 

The  derivation  of  equation  1  to  3  assume  an  average  element  tem¬ 
perature  along  the  perimeter  of  the  tube  for  element  surface  radiation  into 
space.  More  rigorous  analysis  has  shown  that  for  the  specified  element, 
the  computed  numerical  values  can  be  regarded  as  2%  higher  than  the  true 
values.  It  was  also  found  that  inclusion  of  the  internal  radiation  between 
the  inside  surface  of  the  element  would  decrease  the  computed  values  in 
the  same  order.  Both  assumptions  suggest,  therefore,  the  introduction  of 
a  correction  factor  of  0.96  for  average  temperature. 

4.2  Neglecting  the  Overlap 

A  closed  tube  of  equal  wall  thickness  has  been  assumed  for  this 
analysis,  (i.e.  -  equal  conductivity  along  the  cross  section.)  The  over¬ 
lapped  tube  however  has  cross  sections  of  varying  thickness  and  undeter¬ 
mined  thermal  content.  Difficulties  arise  when  the  overlap  is  included  in 
the  equations,  due  to  the  fact  that  the  overlapped  sections  can  have  multiple 
positions  with  respect  to  the  radiation  source.  They  will  change  their 
position  continuously,  due  to  twisting  effect  of  very  long  elements. 

The  evaluation  of  four  distinct  configurations  suggest  that  a  180® 
overlap  will  decrease  the  temperature  gradient  in  the  order  of  25%, 
calling  for  a  further  correction  factor  of  0.75. 
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Suggestion  »f  a  Tolerance  Band  for  the  Results 


..  are  included,  all  the  computed  numerical 

values  of  this  analysis .  would ^be  about ^ 

emissivity  ratios  ,  as  well  *  be  larger  or  smaller  than  assumed, 

the  thermal  expansion  coefficient,  m  y  B  ithin  a  +  30%  band 

In  view  of  this  it  is  felt  that  the  actual  profiles  -educed 

*bh°“  'TTTt  '3T:  <0 are  more  accurately  determined. 
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PART  II 


MATHEMATICAL,  ANALYSIS 
1.0  INTRODUCTION 

The  purpose  of  this  part  of  the  report  is  to  show  the  mathematical 
development  of  the  formulas  and  equations  needed  to  calculate  distorted 
shapes  of  long  thin-walled  tubes  subjected  to  heat  radiation  in  space.  For 
the  most  part,  the  equations  are  developed  in  general  form,  but  in  certain 
cases  where  complexities  warrant  it,  numerical  values  are  used  to  arrive 
at  working  equations.  In  these  cases  the  equations  are  made  to  apply  to  a 
0.5  inch  diameter  beryllium  copper  tube,  0.002  inches  thick,  subjected 
to  radiation  in  space  just  outside  the  earth’s  atmosphere.  Thus  the  prin¬ 
ciple  sources  of  radiation  considered  are  the  sun,  which  delivers  heat 
at  an  intensity  of  450  BTU/HR-FT2,  and  the  earth,  which  delivers  its  own 
heat  at  an  intensity  of  70  BTU/HR-FT2  and  reflected  sun  heat  at  an  intensity 
of  160  BTU/HR-FT2. 

The  method  of  obtaining  the  required  equations  is  to  develope  temper¬ 
ature  profiles  on  the  tube  cross-section  by  consideration  of  the  various 
radiation  and  conduction  heat  flows.  The  temperature  profiles  are  then  used 
to  determine  the  thermal  distortion  of  the  tube. 

Since  the  elements  are  made  of  flat  metal  strips  formed  into  circular 
cross-section  tubes  with  approximately  180  degrees  of  overlap,  it  is 
assumed  that  there  is  always  some  area  of  contact  in  the  overlapped 
section.  Because  of  the  forming  method  used  in  producing  the  tubes,  the 
most  probable  place  for  contact  to  occur  is  at  the  outer  free  edge  of  the 
formed  strip.  This  type  of  contact  makes  the  overlapped  tube  similar  to 
a  seamless  tube  from  a  heat  transfer  point  of  view,  and  for  this  reason 
the  temperature  profile  development  in  section  2.0  deals  only  with  seamless 
tubes.  The  effect  of  increased  area  of  contact  in  the  overlapped  section  on 
temperature  profiles  is  discussed  in  section  5.0. 

Equations  for  the  distorted  shape  of  tubular  elements  are  developed 
in  section  3.0.  These  equations  are  of  interest  mainly  because  they  show 
clearly  what  factors  will  increase  or  decrease  element  deflections .  Actual 
magnitudes  of  deflections  obtained  from  these  equations  are  only  approxi 
mate  because  of  the  assumptions  used  in  the  mathematical  development. 
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2.0  DERIVATION  OF  TEMPERATURE  PROFILES  FOR  THIN-WALLED 
CYLINDRICAL  TUBES 

2.1  Heat  Balance 


The  temperature  profile  of  a  thin -walled  cylindrical  tube  can  be 
determined  by  equating  the  heat  absorbed  to  the  heat  rejected  for  a  dif¬ 
ferential  element  of  the  tube.  Since  the  thickness  t  is  much  smaller  than 
the  radius  r,  (less  than  0.01  r)  temperature  gradients  in  the  radial 
direction  will  be  negligible  compared  to  temperature  gradients  around  the 
perimeter,  and  the  differential  element  shown  in  Figure  1  (a)  can  be 
changed  to  a  rectangular  element  of  length  dx  =  rd/3  ,  as  shown  in 
Figure  1  (b).  The  heat  flows  indicated  in  Figure  1  (b)  are  defined  as 
follows: 


qs  =  Radiant  heat  from  sun  or  other  source 
absorbed  at  the  surface  of  the  tube. 

qr  =  Heat  radiated  to  space. 

qc  =  Heat  flow  around  tube  due  to  conduction. 

q^  =  Heat  radiation  across  interior  of  tube. 

Heat  flow  along  the  length  of  the  tube  is  assumed  to  be  zero,  thus 
allowing  all  further  calculations  to  be  based  on  a  unit  length  of  tube.  This 
assumption  is  valid  as  long  as  the  heat  source  remains  constant  over  long 
lengths  of  the  tube. 
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Equating  heat  into  heat  out  yields 

-Ur*~  fr+ft~  fr*  (1> 

The  temperature  profile  can  be  obtained  from  equation  (1)  by  sub¬ 
stituting  for  the  indicated  heat  flows. 

2.2  Heat  Flow  From  Source  Absorbed  at  Surface 

If  the  source  of  radiant  heat  is  the  sun,  the  heat  flow  is: 

)  <2> 

where  Jfl  is  the  radiation  intensity  and  e*a  is  the  solar  absorptivity  of  the 
material. 

If  the  source  of  radiant  heat  is  the  earth,  the  heal  flow  is: 

f*  ~(jh«s  +■  ces(?j  (3) 

where  jRis  the  intensity  of  sunlight  reflected  by  the  earth,  JK>s  the 
intensity  of  infra-red  earth  radiation  and  £  is  the  emmissivity  of  the  tube 
material.  It  is  assumed  the  absorptivity  of  any  material  for  earth  reflected 
sunlight  is  the  same  as  that  for  direct  sunlight. 

2.3  Heat  Radiated  to  Space 

If  a  surface  radiates  to  black  body  (space)  surroundings,  the  heat 
flow  is: 


j.r  =  reT *4*.  (4) 
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2*4  Heat  Flow  Due  to  Conduction 


from  which: 


(5) 


2,5  Heat  Radiated  Across  Interior  Of  Tube 

i  An-Tder  thC  interior  surface  of  a  tube  divided  into  n  elements  of 

length  Ax  as  shown  in  Figure  2.  The  amount  of  radiation  leaving 
element  1  and  arriving  at  any  element  p  is: 

&  =  j r-W/ 

/  '  r  I  TT  y* 


Figure  2 


where  Jj  is  the  surface  radiosity  of  element  1.  Since  the  elements  lie 
a  circle. /3,  =  /3/>  ,  and  .  Therefore, 


on 


Similarly: 


Ihus  the  net  heat  flow  from  1  to  p  is: 
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And  the  total  heat  flow  from  1  to  all  the  other  elements  is: 

/>  -[(’-tt-iJ-JM  <?> 

Another  expression  for  the  heat  leaving  element  1  can  be  derived 
from  the  fact  that  it  must  be  equal  to  the  negative  sum  of  the  heat  flows 
leaving  the  other  elements.  I.E.: 

$r-o 

Using  the  general  expression  for  heat  leaving  the  surface  of  element  p 


A-##!,-*) 

where  E^  is  the  black  body  emissive  power  of  the  element  and  S'  is 
the  reflectivity,  equation  (8)  becomes 

Combining  equations  (7),  (9),  (10),  and  realizing  that  Efa  = 
gives  the  following  result:  P 

f. 


(10) 

<r^Tp4- 


(ii) 
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If  the  number  of  elements  is  large  such  that  approaches  zero,  then 

(n-jAx  ~  =  ird 

and 

re  | 

vwhere  Tm  18  the  mean  heat  transfer  temperature  of  the  tube.  Thus  the 
required  expression  for  radiation  across  the  interior  of  the  tube  becomes: 


2-6  General  Heat  Flow  Equation 


Substitution  for  the  various  heat  flows  in  equation  (1)  give 
general  heat  balance  in  terms  of  temperatures ,  from  which  the 
profile  may  be  obtained.  Considering  only  sun  radiation  for  the 
equation  (1)  becomes: 


s  the 

temperature 
moment , 


<«(*-)/* 

from  which 


800°  R  and  that  mean  tube  temperature  will  be  around 

T  -  T  *  ?,no‘  aC‘Ual  *«“!»•«*»•  wlU  have  a  maximum  range  of 

from  Post  .  '  ?*“  be8®en  that  th«  function  of  (T/T  )  *  will  range 
from  h 05  to  1  while  cos  (x/r)  ranges  from  1  to  0.  Therefore  to  simplify 
equation  (13)  so  that  an  analytic  solution  may  be  obtained,  it  will  be 

«■«*  *»  •  ■  h.« 
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Stated  in  another  way,  equation  (14)  neglects  the  heat  flow  across  the 
interior  of  the  tube,  and  assumes  that  the  heat  radiated  to  space  from  the 
tube  surface  does  not  vary  from  point  to  point  around  the  perimeter. 

2.7  Temperature  Profile  with  Sun  Radiation  Only 

If  the  sun  is  the  only  source  of  radiation  impinging  on  the  tube  sur¬ 
face,  the  profile  may  be  obtained  from  two  differential  equations.  One 
equation  gives  the  profile  in  the  first  quadrant,  the  other  gives  the  pro¬ 
file  in  the  second  quadrant.  Figure  3  shows  that  by  symmetry  the  profiles 
in  the  third  and  fourth  quadrants  are  identical  to  those  in  the  second  and 
first  quadrants  respectively. 


Figure  3 


The  two  equations  may  be  written  directly  from  equation  (14) 
realizing  that  the  source  is  zero  in  the  second  quadrant. 


[¥<*<Tr] 
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“1 


t 


Equations  (15)  can  be  integrated  directly  to  give: 

#  =  *Tl 

tl  -  f  q,  (?)  *  cj 


,  4Th!,C°”Bta"t8  °f  irlteS™tion  in  equations  (16)  can  be  solved  for  by 
applying  tbe  boundary  conditions  of  continuous  heat  flow  and  temperature 

^tlcaUyt  ‘  2Cr°  fl°W  at  *  =  °  and*  =  Stated  math- 


*  ’Htfal'r)  a  O 
V(*-¥)  -  T,(*  ■  rj 


from  which: 


V  =  c]  1,71 


The  unknown  constant  C  can  he  seen  by  inspection  to  be  the  lowest 

m^nTlt?  °n  f  CrOSS  8eCtlon-  11  can  be  evaluated  in  terms  of  the 
mean  heat  transfer  temperature  by  means  of  the  following  expression- 
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Because  the  variation  of  Tj  and  T  is  small  compared  to  the  value  of 
T^  ,  this  can  be  approximated  by: 


Using  equation  (18)  to  evaluate  C  ,  and  substituting  C  back  into 
equations  (17)  gives  the  desired  temperature  profile.  A  further  substi¬ 
tution  can  be  made  from  an  overall  heat  balance  of  heat  absorbed  to  heat 
radiated: 


Js*s  -  irr&T? 

The  equations  for  temperature  therefore  become: 

7 r  *  ^  *  'kijjkfrT* cos(f)  -  0  4-sJ 

Tt  w  ^  +i£lifc(r-t)*-A4s]  (19) 

Inspection  shows  that  the  maximum  temperature  always  occurs  at 
x  =  0,  and  the  minimum  at  x  =  TTr.  The  equation  for  maximum  temperature 
difference  iB  therefore: 


2.8 


Temperature  Profile 


with  Sun  and  Earth  Radiation 


Directly  Opposed 


If  sun  and  earth  radiation  are  directly  opposed,  the  problem  may  be 
solved  by  two  differential  equations,  since  the  profile  will  still  be  sym¬ 
metrical  as  shown  by  Figure  4. 
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The  two  equations  may  be  written  directly  from  equation  ( 14)  as 
follows: 


0  =  -Tt-^&-C0S(r)  [°<*<¥]  us) 

0  *  [¥<*■<”] 


Replacing  Jg  o(  g  by  Sj  ,  and  (  Jr0<s  +  J£  £ 


R°VS  T  JEC  '  by  Sz‘ 

the  temperature  profile  may  be  obtained  by  direct  integration  of  equation. 
(15),  the  boundary  conditions  being  exactly  the  same  as  with  sun  radiation 
only.  The  desired  profile  is  therefore: 


Tt  "  '^[ilrfr)  f  ~  °-+f]  ( i6) 

5  * 


It  is  not  obvious  from  equations  (16)  where  the  maximum  and  min- 
imum  temperatures  occur  in  the  profile,  the  positions  depending  to  a  large 
extend  on  the  relative  value,  of  S.  and  S  .  To  determine  the  maxima  “ 
and  minima,  equations  (16)  must  te  differentiated  and  set  equal  to  zero- 
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*  0  =  *> 

-£T  =  0^~4f(t~  ty+rfs^y^fr)  b> 

from  which: 

■fs  [“IT  SiJ$y($)]~  /  ^**5 

-t  =  >«(})/£)]-  /  = 

Using  the  values  of  Jg,  J  and  J^,  which  occur  at  the  outer  edges 
of  the  earth's  atmosphere,  valifts  of  o(.3  / €  can  be  obtained  from 
equations  a)  and  b)  for  various  values  of  (x/r). 


-f4-  «  0-087 

a) 

M-f-  j 

H 

f1  *  0-24-Z 

a)  b) 

i &■)-■*  ■, 

■f-  *  TT-I 

-f4-  “  A  39 

b) 

The  significance  of  these  results  is  that  they  give  the  range  of  values 
of  Ofj/jf  for  which  the  minimum  temperature  lies  in  a  certain  quadrant. 
From  the  boundary  conditions  of  the  problem,  it  is  known  that  the  maximum 
temperature  must  be  at  either  x  =  0  or  x  =  *ff  r.  The  result  of  equation  b) 
at  x  =  IT r  shows  that  the  minimum  temperature  will  be  at  x  =  T r  for 
values  of  0<s  /£  ^  1.39,  which  is  the  range  that  covers  most  polished 
metal  surfaces.  It  can  be  seen,  therefore,  that  the  maximum  temperature 
difference  across  a  tube  for  directly  opposed  sun  and  earth  radiation  will 


be:  - 


i 


DHC-SP-TN  164 


=  Tj'fr**))  -  7^(jc*Tr) 

AX*,  =  -&«,+££%  “7| 


2.9 


Temperature  Profile  With  Sun  and 
Opposed 


Earth  Radiation  Not  Directly 


If  the  tube  is  oriented  such  that  the  sun  and  earth  radiation  directions 
have  an  angle  i  between,  the  problem  of  deriving  the  temperature  profile 
becomes  more  complicated.  Four  differential  equations  are  required, 
corresponding  to  the  four  segments  shown  in  Figure  5. 


FLmre  J 


*5 


The  four  equations  may  be  written  directly  from  equation  (14)  as 
follows: 


/T  m  nrJ 

A?  kt 


EHte-f]  08) 


0 
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1 


=  [?<£<r+f] 

-0  -  -fF  [*X*<f] 


Equations  (18)  may  be  integrated  directly  to  give  the  temperature 
profile,  the  boundary  conditions  being  continuity  of  temperature  and 
heat  flow  between  the  segments . 

X  -  WrM' 

X  ~  ir^0)l-X(r)  ~si(h  $r) * c] 


To  obtain  the  maximum  temperature  difference  across  the  cross 
section  of  the  tube,  the  positions  of  maximum  and  minimum  temperature 
must  first  be  established  by  differentiating  equations  (19)  and  setting  the 
derivatives  equal  to  zero. 


4?  *  0  *  *T0)-S,  0  f)  -  S,  (4) 
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Sr'0' 

it*  iJ=  ~l¥A(r)  -  s,  -  S1(l  t  ^ 


V«,  *  S2-  These  J ^yZ eon J "**«■  Sl«>  *  S, 
i*  greater  than  3,  which  is  the  case  that  •Mv*'1  Con8tant  “the  ratio  Kj/i 
be  seen  from  Figure  5  and £ b*  Croidered  W*-  X*  “»  £»•• 
temperature  must  occur  in  segment  1  or  ?  "j"18,  ' 8  that  t*le  maximum 

3  or  4.  Figure  6  shows  the  functions  (x/r)  ’  *"  I"1/11."*"*1  in  8e8ment 


By  aubfltituting  values  of  (x/r)  .  .  . 

(20),  the  maximum  temperature  difference  aero..  0?in/  1“*  ln*°  equationa 
for  any  value  of  f  .  in  .  Ce  acr08s  “•  (ube  can  be  obtained 

difference  would  be  of  the  form:  ’  ®Xpre88lon  for  thl»  temperature 


(21) 
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where  £  indicates  an  unknown  functional  relationship.  Taking  Sj  and  S2 
to  be  constant  in  this  analysis  ,  A  Tmax<  can  be  plotted  as  a  function  of 
f  .  Figure  7  shows  that  the  functional  relationship  approaches  a  cosine 
squared  curve,  such  that  temperature  difference  across  the  tube  can  be 
approximated  by  the  formula: 


“  w(si  +  4 

[°«<f]  (22) 

AX**  *  1 

For  an  element  on  a  satellite  orbiting  the  earth,  the  assumption 
that  is  a  constant  is  not  correct,  since  it  will  change  with  distance 
from  the  earth's  surface  and  with  position  relative  to  the  earth's  bright 
side.  However,  numerical  computations  for  differing  values  of  show 
that  equation  (22)  is  reasonably  accurate  regardless  of  the  value  of  S2, 
in  the  region  52*90°. 


3.0  DERIVATION  OF  THERMAL  DISTORTION 

3. 1  Thermal  Distortion  Due  to  Sun  and  Directly  Opposed  Sun  and  Earth 
Radiation  "  ~~™~~ 

It  has  been  shown  that  a  thin-walled  tube  in  space  will  develops  a 
temperature  difference  across  its  cross-section  due  to  radiation  incident 
on  its  surface.  Because  of  this  temperature  difference,  longitudinal  strips 
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of  the  tub*  will  expand  to  different  lengths,  causing  distortion  of  the  tube. 
Consider  a  short  length  of  such  a  tube,  as  shown  in  Figure  8.  If  the  tem¬ 
perature  gradient  across  the  tube  is  linear  with  respect  to  a  diameter 
passing  through  the  points  of  maximum  and  minimum  temperature,  the 
distortion  will  be  uniform  bending  about  a  radius  of  curvature  R  ,  and 
there  will  be  no  distortion  of  the  cross-section.  c 


Figure  8 


From  the  geometry  of  Figure  8 


From  which 


where  e  is  the  coefficient  of  thermal  expansion  of  the  tube  material  and 
At  ii  the  temperature  difference  across  the  diameter. 

Consider  now  a  long  tube  fixed  at  one  end  with  radiation  striking  it 
at  an  angle  to  the  fixed  end,  as  shown  in  Figure  9. 
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I 


Figure  9 

If  heat  flow  along  the  length  of  the  tube  is  neglected,  the  temperature 
difference  across  the  tube  at  any  angle  6  is: 

AT  »  AT\„  sin  6  (24) 


where  Z\T_  is  the  temperature  difference  at  the  fixed  end. 
max. 

Considering  an  element  of  length  dl  at  any  angle  6 


-7 JU 


(25) 


(26) 
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where  N  =  e  A  T  U. 
dition  that  at  y  =  ; 


Integrating 


equation  (26)  and  applying  the  con- 


-Ny,  (27] 

Taking  the  contangent  of  equation  (27)  gives: 

“  cotfa  -  V y) 

»  C0t(&t-Ny)Ai  (28) 


ate  «  4s- 
d*. 


Integrating  equation  (2  8)  and  applying  the  condition  that  at  y  =  o,  x  =  o; 


* 


sLn6a 


(29) 


Equation  29)  is  the  equation  for  the  shape  a  tube  will  take  .i  radiation 
causes  a  linear  temperature  gradient  along  a  diameter.  Since  0.  can 
have  any  value  from  0  to  180  degrees,  it  is  useful  to  take  9  =  V /2 
in  which  case  equation  (29)  becomes  1 


“  7v"  hi  cos, 


h) 


(30) 


Equation  (30)  can  then  be  plotted  as  shown  in  Figure  5,  Part  I  of  this 

4x18  P°aitive  in  either  Erection  from  zero,  and  the 
*ng1;*/  *”<*icated  h*  aPProPriat8  Points  on  the  curve.  In  this  way, 
e  factor  N  is  calculated  from  the  value  of  ^  T  obtained  bv  con 
eidering  the  radiation  to  be  at  right  angle,  tofhe  tibe,  * 
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Distance  along  the  curve  in  unit  lengths  of  tube  is  also  plotted  in 
Figure  5,  Part  I.  This  is  obtained  by  combining  equations  (24)  and  (25) 
in  a  slightly  different  manner  to  give 


eU---j!j-CSC&4£  (31) 


Integrating  equation  (31)  and  applying  the  condition  (from  Figure  5,  Part  I) 
that  at  1  =  0,  6  =  Til 


L  "  ~ir  l*  tan -fr- 


For  short  tubes,  Figure  5,  Part  I,  shows  that  deflections  are  small  com¬ 
pared  to  the  tube  length.  Thus  the  ordinate  y  in  equation  (30)  may  be 
replaced  by  distance  1  along  the  tube,  giving  the  deflection  x  as 


X. 


(33) 


In  particular  if  N1  ^  0.3,  equation  (33)  simplifies  to 


(34) 


with  a  deviation  of  less  than  2%. 

Equation  (34)  or  Figure  5,  Part  I,  may  be  used  to  give  the  distortion  shape 
of  a  tube  with  a  temperature  difference  across  its  cross-section  induced 
either  by  sun  radiation  or  by  directly  opposed  earth  and  sun  radiation. 
However,  it  must  be  remembered  that  the  shape  is  only  approximate, 
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since  actual  temperature  gradients  with  respect  to  a  diameter  are  not 
linear  as  was  first  stipulated. 

♦  /18U!e  11i8how8  8Qme  Epical  temperature  profiles  plotted  against 

tube  diameter.  It  is  expected  that  the  small  deviations  from  a  linear 

fhatleni  WllluCaUSe,dlSt°rtl0n  °f  the  cr°8s -section,  thus  changing  somewhat 
requir'e^the  8haP<!  ^  ^  T°  CalcuUt*  the  final  Ending  shape  would 


T~  TautL 

&T~rr)4x 


w  ot  0+  H  0'$  ID 


*  **  6*4  U  6*6  hO 
a 


Figure  1 1 


th  n  3n  f°  e"gthwl8e  a"d  cros8"8ectional  strain  energies  to  determine 
the  shape  of  least  energy.  However,  this  final  shape  would  not  be  markedly 
different  from  that  already  derived.  Therefore  the  task  of  deriving  it  will 
not  be  undertaken  in  this  report.  8 

3 ' 2  Iilg.r™al  Distortion  Due  to  Sun  and  Earth  Radiation  Not  Directly 
Opposed  - - - - - - 7 


‘he  sun  and  earth  radiations  are  not  directly  opposed,  the  problem 
•  f  fu  termmln8  the  distortion  shape  becomes  more  complicated.  Not  only 
is  the  temperature  gradient  non-linear  along  a  diameter,  the  points  of 
maximum  and  minimum  temperature  do  not  lie  on  a  diameter.  However 

^  8i.mplified  BOm®what  by  assuming  that  the  temperature 
gradient  due  to  a  single  source  ic  linear  along  a  diameter.  Since  temperature 
profiles  are  solutions  of  linear  differential  equations,  the  profile  due  to 
o  sources  with  an  angle  S  between  is  the  sum  of  two  profiles  due  to 
the  sources  acting  separately.  Also,  since  the  separate  profiles  are 
assumed  linear  the  temperature  difference  across  the  tube  in  the  direction 
of  the  sun  radiation  is  proportional  to  3!  +  S2  cos  5  (see  figure  12)  and 

the  temperature  difference  at  right  angles  to  the  direction  of  sun  radiation 
is  proportional  to  S2  sin  S  .  rauiauon 
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In  this  manner  any  problem  involving  two  radiation  sources  can  be 
transformed  into  a  problem  of  two  different  sources  acting  at  right  angles 
to  each  other . 


Consider  now  a  differential  length  dl  of  the  tube  in  space,  as  shown 
in  Figure  12.  The  effective  source  S3  will  cause  bending  of  the  tube  about  a 


Figure  IJ 


radius  of  curvature  RC3  in  the  direction  of  S3  sin  $  3 .  Similarly  S4  will 
cause  bending  about  a  radius  RC4  in  the  direction  of  S4  sin  6  4.  Using 
equations  previously  developed: 


J-  si„ft 


(35) 
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( 


cJiS 


where 


Integrating  equations  (35)  and  applying  the  conditions  that  at  y  =  z  =  o, 

$3  =  ©4  =  1T/2; 


(36) 


where 


A/ - C  ^3  max 

Ni  ~  cL 


Equations  (35)  may  be  integrated  in  terms  of  the  length  1,  and  then 
solved  for  in  terms  of  $ j  and  6^  to  give: 


ft,  ■  2  arctan Je.  (37) 

&4  -  2.  arctan 

Equations  (36)  and  (37)  may  be  used  to  solve  for  the  components  of  deflection 
(y  and  z)  at  any  distance  1  along  a  tube  which  is  being  deformed  by  two 
sources  of  radiation  with  an  angle  S  between,  each  initially  at  right 
angles  to  the  tube.  The  shape  of  the  deformed  tube  is  shown  in  Figure  14, 


i 
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Figure  14 


Again  it  must  be  pointed  out  that  the  equations  for  y  and  z  give 
approximate  deflections  because  of  the  assumptions  stipulated  in  the 
derivation. 


4.0  DERIVATION  OF  ELEMENT  RESPONSE  TIME 

An  estimate  of  the  length  cf  time  required  for  an  element  to  build  up 
a  temperature  grad  ent  across  its  cross-section  can  be  obtained  by  inves¬ 
tigating  the  response  of  the  mean  tube  temperature  to  sudden  exposure  of 
the  tube  to  sun  radiation.  Tnis  assumes  that  the  temperature  gradients  cal¬ 
culated  in  Section  2.0  will  be  realized  when  the  element  is  in  thermal 
equilibrium. 

Consider  an  element  initially  at  a  uniform  temperature  T  .  If  it 
is  suddenly  placed  in  sunlight,  a  balance  of  the  heat  flows  can  bTSnade  as 
follows: 

Heat  Radiation  In  —  Heat  radiated  out  =  Heat  absorbed . 
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1 


I 


Substituting  for  the  heat  flows  gives  the  following: 

dJsas  -  TcLreTS  =  tyofvrcLt^ 


(38) 


W^reisCHmehe  8PBCifiC  °f  tHe  materU1'  *P  is  the  overlap  factor  and 


Equation  (38)  can  be  re-arranged  to  allow  numerical  integration.  Sub¬ 
stituting  for  the  sun  radiation  in  terms  of  the  final  equilibrium  mean 
temperature  Tm^  gives  the  required  expression: 


(39) 


Fi<m  1  pTerTatUr  Pl0tted  a8  a  fUnCtion  °f  time  ^  i8  ^own  in 

g  re  10,  Part  I,  of  this  report,  for  a  temperature  increase  (T  -  T 
of  285°F.  mf 


mr 


If  the  step  increase  in  temperature  is  small  such  that  T  =  T 

«t*  i  a  _  m  mo 


.  - - - C  AO  BUiail  SUCf 

where  A  Tm  is  very  small  compared  to  T 
equation  (39)  can  be  simplified  to: 


then 


Since  A^m  is  small, 


and  equation  (40)  reduces  to: 
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The  solution  to  equation  (41)  is: 


where 


r  - 

*  4rcT? 


(42) 


and  is  referred  to  as  the  time  constant  of  the  element. 


A.  ti  ^  °VHe  !lmC  constant  C*  i8  of  interest  since  it  represents 

Sure  »earqU1d  l  eiCment  t0  rCaCh  63'2%  °f  itB  temper- 

^derL  ^  '  "  “  *  8tabilitV  analysis  of  an  antenna  which 

undergoes  a  small  change  in  orientation  witn  respent  to  the  radiation  heat 


5.  0  DISCUSSION  OF  ASSUMPTIONS 

5' 1  Interior  Radiation  and  Surface  Radiation  to  Space 


Two  assumptions  made  in  the  derivations  of  the  temperature  profiles 
inte?fr  radiati0n  C0Uld  be  -H'cted  and 
thea‘  *  ‘  •p*S*  at  itB  mean  temperature.  To  determine  the  validity  of 
ese  assumptions  ,  the  heat  flow  equation  will  be  solved  again  for  the 

EquationVl  3 1°*  BUn/adiation  *■*««  all  heat  flow,  into  accost 
Equation  (13)  may  be  rewritten  in  the  following  way: 
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Since  the  difference  between  T  and  Tm  is  very  small,  the  following 
simplification  may  be  used,  with  an  error  of  less  than  2%. 


T'-X.'-lx+fr-V.'fl-X.*  -*ZYr-T-) 


Therefore  equation  (43)  becomes: 


Replacing  the  constants  in  equation  (44)  by  a^  ,  b  and  d,  the  two 
equations  of  heat  flow  for  sun  radiation  only  acting  on  the  tube  may  be 
written  as  follows: 

Jgi i  [r<*<Tr] 

The  solution  to  equations  (45)  is: 


Zafl+a'r*)  sink 


rh&+b&~aa(F) 


Tz  - “2  “(*-?$]  +  a* 

The  maximum  temperature  difterence  across  the  tube  is  therefore: 

-0p. 

■  fkF^s 


(47) 
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where 


It  can  be  seen  that  the  factor  f  is  always  slightly  smaller  than  unity. 
For  the  particular  case  of  a  one-half  inch  beryllium  copper  tube>  0.002 
inches  thick,  f  has  a  value  of  0. 96. 

It  is  therefore  to  be  expected  that  neglecting  interior  radiation  in  the 
tube  and  assuming  radiation  to  space  at  the  mean  temperature  of  the  tube 
will  give  temperature  differences  across  the  tube  which  are  approximately 
4%  too  high. 

5.2  Heat  Flow  Along  the  Tube  Length 

Another  assumption  made  in  the  derivation  of  the  temperature  profile 
was  that  the  heat  flow  along  the  length  of  the  tube  is  zero.  However,  it  can 
be  seen  that  because  of  the  distorted  tube  Bhape ,  the  mean  temperature 
varies  along  the  tube  length  according  to  the  formula: 

T?-T+*U,9  '«> 


where  Tj_ is  the  mean  temperature  at  the  fixed  end,  and  9  is  the  angle  at 
any  point  on  the  tube,  as  shown  in  Figure  9. 

To  obtain  the  average  heat  flux,  equation  (48)  must  be  differentiated 
with  respect  to  length  1  and  multiplied  by  the  conductivity  k  : 
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( 


Substituting  for  d0/dl  from  equation  (31)  gives: 

e0s&  (50) 


Equation  (50)  can  be  differentiated  with  respect  to 6  and  set  equal  to 
aero  to  find  the  point  of  maximum  heat  flux  along  the  length.  If  this  is 
done,  it  is  found  that  ^max  equals  26.6  degrees,  and  the  expression  for 
maximum  heat  flux  becomes: 

(-^Lr°-mkNT"< 


The  average  heat  flux  around  the  perimeter  of  the  tube  at  6  equals 
26.6  degrees  is: 


(52) 


Taking  a  ratio  of  equation  (51)  to  (52),  and  substituting  for  N,  gives,  for 
the  beryllium  copper  tube; 

flux  along  tube  length  =  0.64  eT  'S'  .005 

flux  around  tube  perimeter  m° 

This  result  indicates  that  the  heat  flux  along  the  length  of  the  tube  is 
less  than  one-half  percent  of  the  average  heat  flux  around  the  perimeter  of 
the  tube,  and  will  therefore  have  a  negligible  effect  on  the  temperature 
distributions  calculated  in  section  2 . 

5.3  Overlapped  Tubes 

In  deriving  temperature  profiles  in  Section  2 ,  it  was  assumed  that 
the  tube  was  seamless  and  of  constant  thickness.  If  the  tube  is  overlapped 
180  degrees,  this  assumption  implies  that  there  is  enough  contact  in  the 
overlapped  section  to  give  a  conduction  heat  flow  area  equivalent  to  the 
constant  area  of  the  seamless  tube.  The  amount  of  contact  required  in  this 
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ca.e  would  be  email,  as  illustrated  in  Figure  15  a),  If  the «  contact  are. 
were  larger,  as  shown  in  Figure  15  (b),  the  problem  would  be  aindlar  to 
a  seamless  tube  having  one  side  twice  the  thickness  of  the  other.  For  this 
case,  the  temperature  difference 


Figure  15 


across  the  tube  will  be  derived  with  sun  radiation  at  x  =  0,  *7T r/2  and'Tr. 


When  the  sun  radiation  is  at  x  =  0,  the  problem  is  symmetrical  and 
the  two  equations  of  heat  flow  can  be  written  as  follows: 


^<^1  <«> 

•0--^  [?<*<”■] 

Equations  (53)  can  be  integrated  as  in  Section  2.  with  the  same  boundary 
conditions ,  to  give: 
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The  maximum  temperature  difference  across  the  tube  is  therefore: 

(55) 


Similarly  it  can  be  shown  that  when  the  sun  radiation  is  at  x  =  TTr/2, 
the  maximum  temperature  difference  is: 


wM,  -  (-  1561 

When  the  sun  radiation  is  at  x  =  ITr/Z ,  the  problem  is  no  longer 
symmetrical,  and  four  equations  of  heat  flow  are  required  for  the  four 
quadrants  of  the  tuba.  The  heat  flow  equations  are: 


l‘«Y] 

45 -?F  [”<«<¥] 
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The  solutions  to  equations  (57)  are: 

Z '  ■&&*[£&  ^ ' 1 c] 


By  observation,  it  can  be  seen  that  the  maximum  temperature  must 
occur  in  Tj,  and  the  minimum  in  T ^ .  By  differ entation ,  the  maximum 
temperature  is  found  to  occur  at  x/r  =  0.557T  and  the  minimum  at  x/r  - 
1 . 394  IT  .  Using  these  values  of  x/r  ,  the  maximum  temperature  difference 
is: 


r1 

1 tfxwrrJ-’tfr/Wr)*  ATm  =  0>  7 (59) 

Thus  it  can  be  seen  that  complete  contact  between  layers  in  an  overlapped 
tube  decreases  the  maximum  temperature  difference  across  the  tube  by 
approximately  25%,  for  sun  radiation  only. 

Further  cases  have  not  been  considered  because  of  the  complexity  of 
the  equations  involved,  but  it  is  felt  that  complete  contact  in  the  overlap 
will  reduce  temperature  differences  by  the  same  order  of  magnitude  for 
all  cases. 
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PART  III 


LABORATORY  TESTS 


1.0  INTRODUCTION 

The  purpoae  of  this  teat  ia  to  meaaure  the  element  tip  deflection  aa 
a  function  of  overlap,  direction  of  incident  radiation  with  reapect  to  overlap, 
and  temperature  profile  on  the  crosa-aection  and  along  the  length.  The 
teat  apparatua  may  be  used  to  inveatigate  the  effect  of  changea  in  tubular 
element  deaign  auch  as  aurface  treatment,  perforation,  etc. ,  giving 
ordera  of  magnitude  of  improvementa  which  may  be  expected  in  a  apace 
environment. 


2 . 0  DESCRIPTION  OF  THE  TEST  APPARATUS 

The  layout  of  the  teat  apparatua  ia  ahown  in  Figure  1  to  6.  A  double 
wall  tube  (1).  9  inchea  inaide  diameter  and  aix  feet  high,  encloaea  the 
upright  tubular  element  (2)  to  be  teated.  Thia  configuration  eliminatea 
the  influence  of  weight  aa  the  expected  tip  deflection  will  be  of  the  order 
of  only  1/2  inch.  The  double-wall  tube  ia  watercooled  and  evacuated  and 
ia  closed  on  the  top  by  a  glass  plate  (3)  with  a  calibrated  grid  (4)  to  measure 
the  tip  deflection.  A  heater  element  (5)  is  mounted  parallel  to  the  tubular 
element  in  a  watercooled  housing  (6).  The  length  of  the  heater  element 
is  almost  5  feet  and  the  capacity  ia  3  KW.  The  housing  of  the  he, iter 
element  ia  of  such  design  that  radiant  energy  is  emitted  normal  to  the 
tubular  element  through  a  3/8  inch  slot  (7).  The  dissipated  heat  can  be 
measured  separately  for  each  cooling  device  by  measuring  the  water  flow 
and  temperature.  The  temperatures  of  the  tubular  element  are  measured 
by  thermocouples  soldered  on  the  bottom,  the  middle  and  the  top  of  the 
element,  both  on  the  bright  and  the  dark  aide. 


3 . 0  TEST  LIMITATIONS 

The  test  apparatus  will  not  allow  determination  of  absorptivity  and 
reflectivity  ratios  for  the  tubular  element  with  respect  to  sun  light.  Thia  is 
because  the  ratios  are  obtained  from  the  amount  of  solar  radiation  absorbed 
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or  reflected  in  unit  time  by  a  unit  area,  integrated  over  the  frequencies 
of  the  total  sun  light  spectrum.  In  the  test  apparatus  the  radiation  ratios 
will  be  different  because  the  band  of  frequencies  of  the  heater  element  is 
different  and  the  frequencies  are  lower.  However,  comparature  tests 
between  materials  and  geometry  are  possible  and  will  provide  useful  data. 
It  can  be  safely  assumed  that  any  improvement  in  tube  behavior  which  is 
noted  in  the  test  apparatus  after  a  modification  in  tube  design,  will  produce 
a  similar  improvement  at  least  in  direction  and  order  of  magnitude  for  the 
element  in  space. 


4.0  TEST  PROCpr’oRE 

Pre.U*r‘.»ary  tests  showed  that  tip  deflections  were  affected  by  the 
thermoc.»-.ple  leads  soldered  on  the  tubular  elements.  For  this  reason  it 
was  d'-ided  that  two  identical  elements  would  be  used  for  test  purposes, 

*-o  allow  temperature  measurement  with  soldered  thermocouples,  the 
>iher  to  allow  tip  deflection  measurements  without  the  thermocouples 
attached. 

It  was  intended  that  no  test  results  would  be  tabulated  until  an 
accurate  heat  balance  could  be  obtained,  however,  this  was  not  possible 
in  the  time  available.  The  best  heat  balance  obtained  showed  an  unaccounted 
for  heat  loss  of  20%. 


5.0  TEST  RESULTS 

5* 1  Temperature  Profile  Along  the  Element  Length 

Temperature  on  the  foot  560  to  570°R 
Temperature  in  the  middle  586  to  591°R 
Temperature  on  the  tip  584  to  589°R 

The  temperatures  are  lower  than  they  would  be  in  space,  which  is  to 
be  expected  due  to  the  different  heat  source.  Although  sufficient  time  was 
allowed  for  the  temperatures  to  settle,  the  gradients  persisted  and  were 
observed  ir.  repeated  tests.  This  suggests  that  the  instrumentation  or  the 
heater  input  circuit  must  be  improved. 
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5.2  Temperature  Gradient  Along  the  Perimeter 

The  maximum  temperature  differences  between  the  bright  and  the 
dark  side  of  the  element  were  found  to  be  as  follows: 

on  the  foot  9°F 

in  the  middle  6.3°F 

on  the  tip  5  •  4°F 

These  values  also  are  lower  than  they  would  be  in  sun  light.  There 
should  be  no  difference  in  the  values  on  the  foot  and  in  the  middle  and  on 
the  tip.  However,  as  the  temperatures  along  the  length  differ  it  is  under¬ 
standable  that  a  similar  variation  takes  place  for  the  temperatures  around 
the  perimeter. 

5.3  Tip  Deflection 

The  deflection  was  measured  to  be  .35  inches.  The  theoretical 
deflection  for  a  tubular  element  with  the  same  temperature  gradient -around 
it's  perimeter  is  .45  inches. 


6.0  DISCUSSION  OF  THE  TEST  RESULTS 

The  limited  test  results  indicate  that  the  accuracy  of  the  temperature 
measurements  must  be  improved  prior  to  further  tests.  It  is  believed 
that  this  can  be  achieved  by  selecting  thermocouple  materials  and  recording 
instruments  that  are  more  sensitive  than  those  presently  in.  use.  The 
temperature  of  the  heater  element  should  be  measured  along;  its  length 
to  help  explain  the  temperature  gradient  measured  along  the  length  of  the 
tubular  element. 
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Part  III 


Figure  4  Reading  Temperatures 


Figure  5.  Measuring  Tip  Deflection 


Figure  6.  View  Through  the  Glassiover 


PART  IV 


CONCLUSIONS  AND  RECOMMENDATIONS 
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PART  IV 


CONCLUSIONS  AND  RECOMMENDATIONS 

1.0  SCOPE  OF  THE  MATHEMATICAL  ANALYSIS 

It  Is  felt  that  the  mathematical  analysis  has  been  developed  sufficiently 
to  provide  a  basis  for  predicting  the  distorted  shapes  of  long  tubular 
elements  when  exposed  to  radiation  in  space.  Although  test  results  to  date 
have  not  completely  verified  the  analysis,  the  equations  can  be  used  to 
predict  the  general  shape  and  the  orders  of  magnitude  of  the  radius  of 
curvature  and  tip  deflection  of  tubular  elements.  The  equations  also 
show  the  influence  of  design  par  ameters  such  that  the  effects  of  design 
modifications  maybe  predicted  quite  accurately. 

2.0  INFLUENCE  OF  DESIGN  PARAMETERS 

If  it  is  desired  to  have  tubular  elements  which  remain  "straight" 
when  exposed  to  radiation  in  space,  then  it  is  necessary  to  choose  design 
parameters  which  will  decrease  as  much  as  possible  the  Profile  Number 

so  that  the  radius  of  curvature  becomes  very  large  according  to  the 
formula: 


and  the  tip  deflection  becomes  small  according  to  the  formula: 


If  the  intensity  of  the  incident  radiation  and  the  required  diameter 
of  the  element  are  known,  then  the  following  parameters  offer  the  possibility 
of  either  increasing  or  decreasing  the  Profile  Number  i 
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d/t  =  ratio  of  diameter  to  thickness 

o(s  =  solar  absorptivity 

e  =  coefficient  of  thermal  expansion 

k  =  thermal  conductivity 

All  four  parameters  are  properties  of  the  element  material,  since 
the  value  of  the  ratio  d/t  depends  on  the  strength  and  fatigure  properties 
of  the  metal  used.  Thus  the  problem  of  reducing  element  distortion  in  a 
space  environment  is  simply  a  problem  of  choosing  materials  with  a 
satisfactory  combination  of  strength  and  thermal  properties. 

To  illustrate  the  range  of  possibilities  some  figures  may  be  given 
for  differing  materials. 

The  material  of  the  specified  tube  is  Beryllium  Copper,  with 
e  =  1.05  x  10"6  and  N  =  3  x  10"^.  A  material  such  aB  Invar  has  a  value 
of  e  =  0.7  x  10"6  which  would  reduce  N  to  0.21  x  10-^,  all  other  factors 
remaining  constant.  However,  the  thermal  conductivity  of  Invar  iB 
smaller,  and  the  absorptivity  could  be  larger,  than  that  for  Berylliuqn 
Copper,  giving  a  net  reduction  of  N  not  quite  so  large  as  first  anticipated. 

The  Profile  Number  N  can  also  be  reduced  by  decreasing  the  solar 
absorptivity  of  the  material.  One  method  available  is  to  silver  coat  the 
tube  surface,  which  would  decrease  from  .45  to  .05  and  would  also 
increase  the  overall  thermal  conductivity  from  44  to  80  for  a  silver  to 
Beryllium  Copper  thickness  ratio  of  1  to  4.  This  would  reduce  the  Profile 
Number  N  from  3  x  10"^  to  .  17  x  10*3  assuming  the  factor  d/t  would  re¬ 
main  constant. 

One  other  method  of  reducing  distortion  not  previously  considered  is 
to  modify  the  design  of  the  tube  by  perforating  the  tube  material,  ^y 
doing  this,  some  of  the  radient  energy  would  pass  through  the  perforations 
on  the  bright  side  of  the  tube  and  be  absorbed  on  the  back  side.  In  this  way 
the  temperature  gradient  around  the  perimeter  would  be  reduced,  thus 
causing  smaller  distortions. 


3.0  RECOMMENDATIONS 

Because  of  the  relatively  large  deflections  calculated  for  the  specified 
Beryllium  Copper  tube,  it  is  felt  that  considerable  effort  should  be  spent 
in  verifying  the  theoretical  calculations  by  extensive  experimental  work, 
and  that  further  tests  be  carried  out  to  determine  the  best  possible  means 
of  overcoming  these  large  deflections. 
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The  verification  of  deflections  calculated  from  the  theoretical 
analysis  can  be  achieved  by  use  of  test  apparatus  such  as  described  in 
Part  ul  of  this  report.  Although  no  attempt  has  been  made  to  duplicate 
sun  radiation  in  this  apparatus,  it  is  sufficient  to  produce  temperature 
gradients  in  the  tube  which  can  be  related  to  the  measured  deflections. 
With  some  modifications  to  the  apparatus  and  test  procedure  it  is  felt 
that  the  emissivity  of  the  tube  material  can  be  determined,  since  the  pro¬ 
cedure  represents  a  controlled  steady -state  heat  transfer  experiment. 

In  order  to  achieve  smaller  deflections  it  is  suggested  that  a  study 
of  materials  be  made  to  determine  which  have  the  best  combination  of 
thermal  and  strength  properties  for  reducing  distortion.  Except  for  solar 
absorptivity,  the  properties  of  most  materials  are  well  known  and  easily 
found.  Thus  the  study  would  require  extensive  experimental  work  in  de¬ 
termining  solar  absorptivities  of  materials  which  might  be  considered 
useful  from  a  consideration  of  the  other  properties  of  interest. 


Combinations  of  materials  to  use  to  best  advantage  the  properties  of 
each  one  should  also  be  considered  and  studied.  One  example  of  this  would 
be  silver  coating  on  a  beryllium  copper  tubular  element.  The  low  solar 
absorptivity  and  high  thermal  conductivity  of  silver  would  combine  well 
with  the  high  strength  properties  ot  beryllium  copper.  Such  a  study  would 
necessarily  require  information  with  regard  to  the  best  means  of  bonding 
the  materials,  optimum  relative  thickness  of  the  materials  and  the  best 
manufacturing  and  finishing  methods  available. 


Methods  for  reducing  distortion  which  primarily  do  not  depend  on 
the  properties  of  the  element  materials  should  also  be  investigated.  One 
such  method  would  involve  perforating  the  tubular  element  so  that  some 
of  the  solar  radiation  might  pass  through  the  hot  side  of  the  tube  and  be 
absorbed  on  the  cold  side,  thus  reducing  temperatures  gradients.  Such  a 
study  would  require  an  extensive  mathematical  analysis  to  determine 
patterns  for  the  perforations  and  also  to  determine  the  ettects  of  the 
perforations  on  radiation  areas,  heat  conduction  paths,  and  tube  strength. 
The  study  would  also  require  experimental  tests  to  verify  any  theoretical 
calculations . 
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APPENDICES 


LIST  OF  SYMBOLS 


REFERENCES 


LIST  OF  SYMBOLS 


units 


A  Satellite  altitude  above  earth  surface 

B  Earth  radiant  energy  factor  (see  3.2,  Part  1) 

F .  Altitude  correction  factor  (see  3.1,  Part  I) 

A 

j  Radiant  energy 

N  Element  profile  number  (see  2.3,  Part  I) 

R  Earth  radius 

R  Radius  of  curvature 

c 

T  Temperature 

Index  S,s  Refers  to  sun  radiation 
Index  E,e  Refers  to  earth  radiation 
Index  R.r  Refers  to  earth  reflection 
c  Specific  heat 

d  Element  diameter 

e  Thermal  expansion  coefficient 

k  Thermal  conductivity 

1  Element  length 

t  Element  material  thickness 

r  Element  radius 


miles 


BTU/hr  .ft2 

ffl 

miles 

feet 

°R 


Btu/lb-°F 
feet 
in /in 

BTU/ft.  °F  hr 

feet 

feet 

feet 


x.y.z  Axes  of  co-ordinate  system 

0(  Absorptivity  ratio 


Ml*  MOT!  .  MM  M. 


jr 

5 

€ 

"( 

9 

<P 

<r 

<* 

q 

Eb 

c 


List  of  Symbols  (cont'd) 
Reflectivity  ratio 

Angle  between  sun  and  earth  radiation 
direction  with  respect  to  the  crossection 
of  the  element 

Emissivity  ratio 

Angle  between  perpendiculars  from  earth 
to  satellite  and  earth  to  sun 

Radiation  incident  angle  -  angle  between 
the  element  axis  and  the  radiation 
direction  seen  from  the  root  of  the 
element. 

Overlap  factor 

Material  density 

Stefan-Boltzmann  constant 

Time 

Element  time  constant 
Heat  flow 

Black  body  emissive  power 
Constant  of  integration 


degrees 

degrees 

degrees 

lb  /ft3 

BTU/ft2hr.(°F)4 

minutes 

minutes 

BTU/hr. 

BTU/hr-ft2 
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